HEPHY-PUB 700/98 
STPHY 30/98 



MULTIPLICITY DEPENDENCE OF LIKE-SIGN AND 
OPPOSITE-SIGN CORRELATIONS IN pp REACTIONS 

B. BUSCHBECK*, H.C. EGGERS 1 " and P. LIPA* 
* Institut fur Hochenergiephysik, Nikolsdorfergasse 18, A-1050 Vienna, Austria 
' Department of Physics, University of Stellenbosch, 7600 Stellenbosch, South Africa 



Abstract 

Discussions about Bose-Einstein correlations between decay products of copro- 
duced W-bosons again raise the question about the behaviour of correlations if 
several strings are produced. This is studied by the multiplicity dependence of 
correlation functions of particle pairs with like-sign and opposite-sign charge in pp 
reactions at *Js = 630 GeV. 
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1 Introduction 

Recently, there has been much discussion regarding the. possibility that Bose- 
Einstein correlations and other interconnection effects □ between decay prod- 
ucts of different strings could affect the measurement of the W mass. Since 
measurements are hampered by low statistics and experimental difficulties, the 
question arises whether and where effects of the superposition of several strings 
can be tested independently. Within the Dual Parton Model, the number of 
strings may be expected to be proportional to multiplicity and thus influence 
the multiplicity dependence of various effects, in particular Bose-Einstein cor- 
relations. While the decrease in observed A as function of nmltiplicity can be 
explained naturally in terras of products of different strings, 013 other features 
such as the radius cannot. □ 

Improvements in experimental analysis techniques, i larger data samples 
and the above theoretical background make it desirable to repeat and extend 
Bose-Einstein analyses with an advanced strategy. In this contribution, we 
investigate like-sign and opposite-sign correlations at different total multiplic- 
ities with the same model-independent strategy and good statistics. The bias 
introduced by selecting events ofj-a given overall multiplicity is eliminated with 
the use of "internal cumulants" . B 



2 Data sample and normalized density correlation functions 

The data sample consists of 600.000 non-single-diffractive pp reactions at s/s — 
630 GeV measured by the UA1 central detector. El As before, only vertex- 
associated charged tracks with transverse momentum px > 0.15 GeV/c, \rj\ < 
3, good measurement quality and fitted length > 30cm have been used. To 
avoid acceptance problems, we restricted the azimuthal angle to 45° < \<j)\ < 
135° ("good azimuth"). 

Fig. 1 shows the normalized density correlation functions for pairs of like- 
sign charge and for opposite-sign charge separately. Restricting the total un- 
corrected charged multiplicity N in |?7| < 3 to the windows 1<A^ < 6 (Fig. 
la) and 28<N < 35 (Fig. lb), one obtains for the corrected particle density 
in the central rapidity region dn/drj = 0.83 ± 0.08 and dn/dij = 5.3 ± 0.17 re- 
spectivelyn. All quantities measured are defined in the notation of correlation 



"Events are selected according to multiplicity in all azimuthal regions, while the subsequent 
analysis is performed for particles in the good azimuth only. This is done because the total 
rather than the good-azimuth multiplicity is physically relevant. Total multiplicity density 
{dn/drj) is estimated as twice the density measured in the good-azimuth region. 
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integrals 0, 

pi{Q) J n d 3 pi d 3 p 2 p 2 (Pi, P2) S[Q- ?(pi, p 2 )] 



r 2 (Q) = 



pi®p\{Q) f n d 3 pi d 3 p 2 Pi(pi) j0i(p 2 ) S[Q- g(pi, p 2 ) 



with p the three-momenta and q = y— (pi — P2) 2 , with p being the correspond- 
ing four-momenta. The integration region £1 is identical with our experimental 
cuts as specified above. All particles have been assumed to be pions. 

In pp reactions and in full phase space the number of positive and nega- 
tive particles are equal; furthermore in the central rapidity region the corre- 
sponding p- functions are also equal. Given the charged multiplicity pi(p) = 
Pl(p) + Px (p)j we hence assume that pf (p) = Pi(p) — |pi(p) and therefore 
also Px®p~\{Q) — jPi®Pi(Q) etc. The like-sign and opposite-sign normalised 
correlation densities hence become, respectively, 

r is (Q) = ptHQ) = pt + (Q) + pr(Q) „ pt + (Q) + pr(Q) 
pt®pf(Q) pt®pt(Q) + Pi®Pi(Q) \pi®pi{Q) 
pf T (Q) _ pi~{Q) + P2 + {Q) ^pt-{Q) + P2 + (Q) 



S {Q) 



pt®pf(Q) pt®pi{Q) + pi®ptiQ) " \pi®px{Q) 



In Fig. 1, both the like-sign and opposite-sign correlation densities show a 
strong dependence on multiplicity. To perform, however, a quantitative anal- 
ysis, one has to get rid of the combinatorial background by calculating the 
cumulants and secondly to correct for the bias introduced by fixing multiplic- 

3 Cumulants for limited multiplicity ranges 

A quantitative study ofJhe bias introduced by fixing multiplicity has been 
performed by Lipa et aLB At fixed total multiplicity N, the standard factorial 
cumulants 

k 2 (Pi, P2 \N) = pa(Pi, Pa \N) ~ Pi(Pi \N) Pl ( Pl \N) (1) 

are nevertheless nonzero, even when particles are completely uncorrelated. 
These purely "external" correlations appear because the cumulant for any 
multinomial distribution is nonzero; for example, in second order, 

^(Pi, P2 \N) = ~Pi(pi \N) Pl (p 2 \N) . (2) 



b The usual Bose-Einstein analysis assumes that r% tends to a constant for large Q, ie. 
r 2 BE 05 — ^) = cons tant. It should be clear from Figure 1 that no such constancy exists for 
limited multiplicity windows. 
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As shown in Ref . El, the "internal cumulants' 



4(Pi,P2 \N) = K 2 ( Pl ,p 2 \N) - K r lt (Pi,P2 \N) (3) 

correct this bias exactly: they are zero whenever the N particles behave multi- 
nomially. Integrating (|l|)-(||) to obtain the correlation integral and normaliz- 
ing, we arrive at the normalized internal cumulants for given fixed N, 

Kl(Q\N)= K ^ Q ) N ) =r 2 (Q\N)- f 1 - — V (4) 

The prescription for calculating the K^(Q \N) is therefore to measure at fixed 
N (in Q) first r 2 (Q \N) = p 2 (Q \N) / pi®px(Q \N) and then to subtract 

i _ 1__ N(N-l) _ J nP2 (Q\N)dQ 



N N 2 S n Pi®Pi{Q\N)dQ' 

To obtain adequate statistics, we have to measure averages over limited mul- 
tiplicity ranges rather than at fixed N . The internaLsecond-order cumulant 
averaged over the multiplicity range [A, B] is given inB by 



— // „ \ En=a p n4(Pi,P2\N) 

^2 (Pl,P 2 ) ~ — g , 



y*— S ■ (5) 



where Pjv is the experimental multiplicity distribution in [A, B] and the bar 
over any quantity S denotes S = J2n=a PnS(N)/ J2n=a Pn- 
Assuming that the shape of pi(p \N) does not vary in [A, B] 

pi(p\N)~(N/N)pl( Pl ), 
and adopting again the correlation integral prescription, we can write 
WUq) = ~ 2l{Q) = MQ) - InP2(Q)dQ 

pi®pi(Q) pi®Pi(Q) f a Pi®Pi{Q) d Q' 



For like-sign measurements, this becomes 



— us = -pj s (Q) JnP2(Q)dQ 

|/Wi(Q) U n Pi^Pi(Q)dQ 
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Figure 1: versus Q for like-sign pairs (full circles) and versus Q for opposite-sign 
pairs (open cirles) at a) dn/drj = 0.83 and b) dn/dn = 5.3. 
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Figure 2: Internal cumulants eqns. (|7j) for is pairs (full circles) and os pairs (open circles) 
for two different multiplicity densities, a) dn/drj = 0.83 and b) dn/drj = 5.3. 
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and similarly for K2 ■ Since all quantities shown here and below are to be 
understood as mean values in multiplicity ranges like in (|5|) - (Q) , we henceforth 
(and in Fig. 1) omit the bar on the symbols. 

Fig. 2 shows both like- and opposite-sign internal cumulants for two selec- 
tions of dn/dn. Three features are immediately apparent: 

I. The cumulants differ from the moments of Fig. 1 by a shift constant 
in Q, but the amount of the shift changes for different cumulants. The 
importance of changing from r2 to K\ lies in the fact that the latter 
demarcate clearly the "line of no correlation" . 

II. Internal cumulants integrate to zero over the entire phase space; hence 
the positive part of K\ at small Q is compensated by a negative part 
at larger Q. Physically, this means that particles like to cluster, so that 
there is a surfeit of pairs at small Q and a dearth of pairs at large Q 
compared to the uncorrelated case. 

III. The dependencies of the like- and opposite-sign cumulants on multiplicity 
are rather similar in that both decrease markedly with dn/drj. This is 
discussed below. 



4 Multiplicity dependence of normalized cumulants 

The similar decrease of K\ and K\ os with dn/drj suggests that both could 
have the same functional dependence on multiplicity density. Under this hy- 
pothesis, 

K^ s {Q,dn/dn) = Y ls (Q)C{dn/ dn,Q) , (8) 
K'° S {Q, dn/dn) = Y os (Q)C(dn/dn,Q) , 

the quotient of the cumulants should be independent of multiplicity 

Ki^dn/dn) _ Y ts (Q) _ f CQnstantm dny 



Kl os (Q, dn/drj) Y os {Q) \ drj 

Figure 3a shows that (^|) holds at least approximately within error bars. (In 
the region where both cumulants are near zero, no meaningful quotients can 
be formed.) 

Having shown that like- and unlike-sign internal cumulants behave approx- 
imately in the same way as functions of dn/dn, we now ask how this dependence 
is structured. A first hypothesis, is that K\ depends inversely on N; in the 
notation of Eq. M) , 
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Ki a (Q\N) = Y a (Q)C(N,Q)=Y a (Q)N- 1 a = £s,os. (10) 
This can be motivated theoretically by 

a) Resonances: If the unnormalized cumulants k\ os and were wholly 
the result of resonance decays and if the number of resonances were 
proportional to the multiplicity A, then k\ oc A. Assuming p±{p \N) oc 
Npi(p) gives pi®pi oc A 2 , and hence after normalization, the resonance- 
inspired guess 

I I CS 1 

R j = ^ _L_ 

pi®pi A 

b) Independent superposition in momentum space of n equal strings 
would also lead to K\ = n ■ K 2 strms /(/9i£g)yOi) oc (1/n) oc (1/N). Devia- 
tions can occur if the strings are unequal. 

Eq. @ implies that K{ a (Q \Ni)/K^ a (Q \N 2 ) = {N 2 /N l ) = (constant in Q) 



for two multiplicities N\ and A^- In Fig. 3b, we show the quotient of cumulants 
for two multiplicities. Surprisingly, we find not one but two constants, one 
for small Q < 0.4 GeV, one for large Q > 2 GeV, where only the latter 
corresponds to the value (N2/N1) expected from (|Tc|), shown as the dotted 
line. At small Q, one must clearly look for othex functional forms for C(N, Q). 
Some phenomenological guesses are as follows. B 

c) In the source picture of Bose-Einstein correlations, the independent 
superposition of sources in configuration space leads to K l ^ BE ~ constant 
in N. We would expect A to be independent of multiplicity, but the 
radius R to increase with N . Hence there should be a change of shape 
of K2 is (Q) but no overall A-dependence. 

d) A mixture of processes, a) with c) could result in a dependence 

K^{Q\N)*a{Q)+ h M (ii) 

However, no comparable source picture is available for the unlike-sign 
case. Also, the quotient K\ jK\ os would not be constant in A, in 
contradiction to the results of Fig. 3a. 

e) Another guess would yield a A -1 dependence for large Q, while the 
small-Q region scales with A~ Q , with a different "constant" Y2, 

Al a = ^^e(Q-2) + ^Me(0.4-Q) + ... a = £s,os, (12) 
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3|) is shown for two selections of dn/dr). b) The ratio of two 

1 has been 

omitted because the denominators are around zero. 



Figure 3: a) The ratio Eq. 

i^| a (Q) at different dn/dr} as indicated on ordinate. The region around 



In order to test the above ideas, we plot the cumulants against {dn/dr]) -1 
as follows. To avoid local statistical fluctuations, the normalized cumulants 
K% (Q) and K2° S (Q) are fitted with suitable functions in restricted Q-ranges 
(not shown). The have also been fitted to an exponential parametriza- 

tion for Q < 1 GeV/c 

K I 2 ls {Q) = a + \e- R Q . (13) 

The multiplicity dependence of the fitted K\ ls and K\ os (circles), as well as 
the corresponding A values (crosses) are plotted in Fig. 4. Fig. 4b shows that, 
as in Fig. 3b, the (l/A f )-dependence is satisfied only for large Q, but not for 
small Q. The a+b/N dependence in Fig. 4a (solid line) provides a possible, but 
hardly unique, explanation. The dashed line corresponding to N~ a behaviour 
( |l2l ) seems somewhat better. However, while this form is compatible with all 



results shown so far, we have no phenomenological justification for it. Eq. (12) 
also omits the region around Q ~ 1 GeV (e.g. the region of p° production) 
where the phase space contributes maximally, but the K\ os decrease rapidly 
with increasing Q and the K\ ls are already small (Fig. 2). It is difficult to 
investigate this important region directly. A X/N dependence (via resonances, 
for example) in the dominant region around 1 GeV/c could presumably cause 
the large-Q region to follow suit via missing pairs. 
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Figure 4: a) Multiplicity dependence of K' 2 ls (Q = 0.11 (full circles) and of K^ os (Q = 0.1) 
(open circles) and A values of exponential fit, Eq. (h3f) (crosses), b) as in a) but for the 
large Q region. For better comparison of the respective dependencies on dn/drj, the absolute 
values have been scaled bv constant factors. The straight lines are linear fits, the dashed 
line corresponds to eqn.(n2h. 



5 Summary 

The multiplicity dependence of like-sign and opposite-sign two-body correla- 
tion functions are studied with the same model-independent strategy. The bias 
introduced by selecting events of a given overall multiplicity is eliminated in 
measuring "internal cumulants" . We observe that 

- the like-sign and opposite-sign functions have a very similar multiplicity 
dependence, which is surprising because of their different shapes and 
assumed origins, 

- there exist two regions, one at small Q, where the multiplicity dependence 
of both is smaller than l/N, and one at large Q( > 1 GeV/c), where the 
functions are negative and follow roughly an 1/7V law. 

Theoretical worlfflll is challenged by these findings. The similar be- 
haviour of Is and os functions at small Q, both decreasing with multiplicity, 
favours suppression of Bose-Einstein correlations between products of different 
strings. Resonance decays (Regge terms) should arguably contribute also to Is- 
functionsBEJ. One could therefore try to explain alternatively the decrease of 
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is functions with N by a mixture of Bose-Einstcin correlations (assumed to be 
constant in N) with resonance production. This leaves unexplained, however, 
the similar behaviour of os and is functions and the nearly TV independent 
shape^j in the small-Q region. 
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